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We evaluate the one-particle and double-occupied Green functions for the Hubbard model at 
half-filling using the moment approach of Nolting Q. Our starting point is a self-energy, £(fc, us), 
which has a single pole, O(fc), with spectral weight, a(k), and quasi-particle lifetime, 7(fc) Q. In our 
approach, £(fc,a>) becomes the central feature of the many-body problem and due to three unkown 
fc-parameters we have to satisfy only the first three sum rules instead of four as in the canonical 
formulation of Nolting Q . This self-energy choice forces our system to be a non-Fermi liquid for any 
value of the interaction, since it does not vanish at zero frequency. The one-particle Green function, 
G(k,to), shows the finger-print of a strongly correlated system, i.e., a double peak structure in the 
■ one-particle spectral density, A(k, u), vs u) for intermediate values of the interaction. Close to the 

Mott Insulator-Transition, A(k,ui), becomes a wide single peak, signaling the absence of quasi- 
particles. Similar behavior is observed for the real and imaginary parts of the self-energy, E(fc, ui). 
The double-occupied Green function, G2 {q, has been obtained from G(k,u) by means of the equa- 
\Q • tion of motion. The relation between G2{q,ui) and the self-energy, Y,(k,ui), is formally established 

' and numerical results for the spectral function of G2(k,uj), x^ 2 \k,uj) = ——lim s ^ + Ira \p2{k, 

are given. Our approach represents the simplest way to include: 1- lifetime effects in the moment 
approach of Nolting, as shown in the paper; 2- Fermi or/and Marginal Fermi liquid features as we 
discuss in the conclusions. 
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I. INTRODUCTION 



After the discovery of the High-X^ materials |j] , the study of correlations has gained interested due to the fact that 
there is the belief that the normal properties of these materials could be explained in the framework of the Hubbard 
model J5|||, since electron correlations are strong, i.e., the on-site electron-electron repulsions U are much larger than 
the energies associated to the hybridization of atomic orbitals belonging to different atoms JtJ . We consider the study 
of correlations in the Hubbard model as a rewarding task since it will shed light on still unsolved points of the novel 
materials. For example, at high temperatures (T c 30 — 130 K) these HTSC cuprates, which are poor conductors, 
become superconductors. This feature is strange indeed because the Coulomb repulsion is strong. Contrary to the 
predictions of the Fermi liquid theory, the resistivity at T > T c and optimum doping is linear in temperature, i.e., 
k T |. This suggests a very strong scattering of elementary excitations. A discussion of the possible breakdown 
of Fermi liquid theory is given in Ref. j^] . In the present work we explore the effects of having a non-Fermi liquid 
behavior into the one-particle and double-occupied dynamical properties. 

We will use the moment approach (or sum rules) of Nolting jlj for the spectral density, A(k, lo). It is well known in 
the literature |i"c| that the moment approach in the spherical approximation - when the narrowing band factor, B(k), 
is not k - dependent - always gives a gap in the density of states (DOS). If the chemical potential happens to be in 
this gap, then we always have an insulator. It has been argued that the way to cure this unrealistic gap is to have a 
better approximation for the narrowing band factor, B(k). 

We have followed a different path which consists in proposing a single pole structure in the self-energy, £(fc,<x>). 
Closing the gap for small and intermediate values of U/W (W is the bandwidth and W = 8t in two dimensions), is 
not the only rationality behind our calculation. It is well documented [p"l| |i~3| that correlations give rise to lifetime 
contributions that conspire against the very definition of quasi-particles at the chemical potential. This leads to an 
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eventual collapse of the Fermi liquid picture of adiabaticity between the electron gas and the electron liquid so well 
described by the phenomenological treatment of Landau. 

This paper is organized as follows. In Section [ll], we present our model Hamiltonian and our self-energy proposal, 
in which we try to justify our bold choice, as one step forward to understand the physics behind the moment approach 
of Nolting. This led us to come up with an Ansatz which is not zero at zero frequency, i.e., our system is not a Fermi 
liquid for any value of the interaction. We could argue that our Ansatz is valid for energy scales not too close to the 
chemical potential. However, our Ansatz has to become exact in the limit of U/W ^> 1. In Section III, we present 
our results which consist of the real and imaginary parts of the self-energy, and spectral functions of the one- and 
two-particle Green functions along the diagonal of the Brillouin zone, for different values of the interaction. We have 
used a system size of 32 x 32. In Section IV we present our conclusions and the future trends. 



II. MODEL HAMILTONIAN AND SELF-ENERGY ANSATZ 



The model we study is the Hubbard Hamiltonian 

H = fr^ct c-! H n^ n^- — acL , (1) 

h3 ia 3° 2 ia 

where ct (cj ) are creation (annihilation) electron operators with spin a. n-? = ct c? . U is the local interaction, 
p the chemical potential and we work in the grand canonical ensemble. We have adopted Einstein convention for 
repeated indices, i.e., for the N s sites i, the z nearest-neighbor sites (n.n.) j and for spin up and down (a = —o = ±1). 
tj* 7 = —t, for n.n. and zero otherwise. 

Let us propose for the self-energy, E(fc, ui), the following single pole Ansatz 



E(£,w)=ptf + ^ =r ; a(fc), 7 (fc)e5ft (2) 

bj — fi(fe) — ij(k) 

With our choice for u>), we have that the real and imaginary parts of the self-energy do satisfy the Kramers- 
Kronig relations |T3[, since it is analytic in one of half of the complex plane. In fact, the physical solution to the 

problem is when a(k) > 0, as it can be checked by finding the roots of G(k,ui) in the complex plane. We postpone 
the discussion of this point for the conclusions. The Ansatz given in Eq. (||) has some similarity with the Hubbard-I 
solution p"5| ]. However, we have neglected any frequency dependence in the damping. Our calculations show that j(k) 
is fc -independent but strongly [/-dependent. 

The validity of Luttinger theorem |l6| has been discussed in Ref. |l7|| . We argue that most likely the Luttinger 
theorem is not going to hold because we have a non-Fermi liquid system. Our choice of pU in Eq. (Q), the Hartree- 
shift, is very convenient since it redefines an effective chemical potential, p e ff — P — pU ■ This effective potential is 
zero at half-filling, p = 1/2, since p = U/2 there. Then, u> = means that we are at the chemical potential. We 
want to explicitly state that our choice (Eq. (||)) has the advantage of requiring only three sum rules to be satisfied, 
instead of four as in the now normal procedure of Nolting jjj which starts from the spectral function, A(k,ui), itself. 
We add that if j(k) = 0, then we go back to Nolting's canonical results. In this case |l8j 

a{k) = p{\ - p)U 2 ; Sl(k) = (1 -p)U + B(k) , B(k) — B + F(k) , (3) 

where B{k) is the narrowing band factor. The fc-independent narrowing band factor, B, is calculated in closed form 
in Ref. It is a self consistent quantity, though. The fc-dependent narrowing band factor, F(k), has been evaluated 
recently by Herrmann and Nolting p9| using a two-pole ansatz with the two poles located at the same energies than 
the poles of the one-particle Green function. This treatment, beyond the spherical treatment of Nolting H, can 
be mapped into the calculations of Kishore and Granato p^] with appropiate identification of our parameters in the 
paramagnetic phase. Their approach gives a Mott Metal-Insulator transition. All this means that the metal-insulator 
transition (MMIT) is embedded into the Hubbard model and it does not require of lifetime effects to accomplish this. 
However, once more, lifetime effects are a natural element of the many-body physics for intermediate and strongly 
correlated electron systems where the concept of quasi-particle docs not apply any longer. We do not pursue anymore 
here the MMIT, but the interested reader is addressed to references p|,p |2l]j22] ] . 

By definition the one-particle Green function, G(k,u>), in terms of ui), is given as 
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where = — 2t(cos(k x ) + cos(k y ) — p + pU. Also, we will require the one-particle spectral density, A{k, uj), which is 
defined as 

A(k,uj) = -— lim ImG(k,uj + iS) (5) 

7T 5^0+ 

Using Eqs. (|-|), we arrive to the following expression for the spectral density 

A(fc,c) = - (6) 



77 ((a; - EjE )(w -%)-«(£)) + 7 2 (fc)(^-e^ 



Using the first three sum rules of Nolting jl| for the the spectral function of Eq. (|^) we obtain the following 
equations 

/+oo _ _ 

A(k,Lu)dLU EE M (k) = 1 
-CO 

+oo 

ujA(k,uj)duj ~ Mi(k) = Ej: 

cu 2 A(k, u)du = M 2 (k) = {er - pU f + 2pU(e £ - pU) + pU 2 , (7) 



oc 
+00 



where the Mj(fc)'s, i = 0, 1, 2, are the first three moments PJ. For example, the first moment (i = 0) is the area below 
the curve of A(k,uj) vs w, the second moment (i = 1) is the center of gravity of the spectral function and the second 
order moment (or third moment, i = 2) is related to the width of the spectral function, A(k, uj). So, damping effects 
are controlled by the second order sum rule. We do not use the fourth moment or sum rule because we have three 
fc-dependent unknown parameters (our way of working is different to one of Nolting since in the latter we have to 
use four moments. The difference lies in the fact that he starts with the one-particle spectral density). We could 
guess that in order to extend the canonical formalism of Nolting to include lifetime effects, starting from his two pole 
ansatz, we should have to postulate the following structure for G(k,uj) 

Gi k,^ = ^ + ^01 ^ , (8 ) 

lo — LOi(k) + i~f(k) uj — uj2(k) + ij(k) 

from where we see that we would need five moments or sum rules because we have five parameters to determine, i.e., 
ai(k), 0Ji(k), 7(fc), with i = 1,2. With the proposal (Eq. (^)) we have only three parameters to calculate. 

We assume that at p — 1/2 the chemical potential, p = U/2. The density of states which results of the two pole 
Ansatz for the one-particle Green function, in the spherical approximation of Nolting |]l|, always has a gap. This 



3 1 which has been critized since many 
where the authors point out to the 



solution (always a gap) is known in the literature as the Hubbard-I solution 
years ago by Laura Roth among others. We call the attention to Ref. 2 
fact that the fc-dependence has to be included in B(k). A recent calculation by Kirchhofer [|26| is performed at the 
mean field level for the fc-dependence of the band narrowing factor, when the two Hubbard bands are separated. In 
a more elaborated calculation based on the Mori's formalism [^7| for the one-particle Green function, Kirchhoffer et 
al [|26U2§| 1 obtain three peaks in the spectral density, A(k, oj), which respects particle-hole symmetry. In the end, they 
get a Mott metal insulator transition, for U/t = 5. Here we are including lifetime effects as a crucial ingredient in 
the formulation beyond a mean-field treatment. Kirchhofer p6| , p8| also considers the presence of antiferromagnetism 
fluctuations in an empirical way. We could extend Kirchhoffer et al's calculations using the numerical values of the 
dynamical spin susceptibility, x(<?j w )j m the spin-fermin model of superconductivity of Pines, Chubukov and others 
p9| . The dynamical spi n s usceptibility has been obtained from nuclear magnetic resonance experiments in the High 



T c cuprates. In Section III, we present our numerical results and their interpretation. 
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III. NUMERICAL RESULTS AND THEIR INTERPRETATION 



In Figs, la, lb and lc we present the spectral density, A(k,co), vs lo along the diagonal of the Brillouin zone 
(k = 27r(n,n)/32) for U/W = 1/2,2/3 and 1, respectively. We are working with a finite system of periodicity of 
32 x 32. For U/W = 1/2 we have a double-peaked structure, with visible lifetime effects (the Dirac delta functions 
of Nolting now have width) . This is a feature of correlated electron systems as it has been discussed in the work of 
Schneider et al jl0| for the case of U < 0. The physics is different but the peak structure is similar. For U/W = 2/3 
we still observe the double peak structure but lifetime effects are stronger. Finally, for U/W = 1, the double peak 
structure is practically washed out. As we see, lifetime effects are very much pronounced for the larger values of 
U/W presented, i.e., for U/W = 1.0. The two peaks of A(k,u>) vs u> are separated approximately by a distance of 

V(^-%-7(fe)) 2 +47(AK-. 

In Figs. 2a, 2b and 2c we show the imaginary part of the self-energy, — Im[T,(k, uj)], vs lo along the diagonal of the 
Brillouin zone for the same values of U/W of Fig. 1. Again we observe that for increasing values of U/W, damping 
effects are stronger in the self-energy, as is the case in the one-particle spectral function (see Figs, la, lb and lc). 
In addition, we do not observe any Fermi liquid dependence (in frequency) of the imaginary part of the self-energy 
around lo = 0. This is due to our choice of our ansatz (Eq. (|2|)). We could include Fermi or Marginal Fermi liquid 



behavior close to the chemical potential as it is suggested in Section IV. However, within the present work, we could 
say that our approximation is valid for frequencies not too close to the chemical potential. Certaintly, for small values 
of U/W, we should have some Fermi liquid behavior (at least in 2-d), like an imaginary self-energy going to zero as 
positive power of lo at the chemical potential |fso| . 

In Fig. 3a, 3b and 3c we report results for the real part of the self-energy, i?e[S(/c, lo)], vs lo along the diagonal 
of the Brillouin zone for the same values of interaction of Fig. 1. For U /W — 1/2, 2/3 we see a more o less regular 
pattern. However, for U/W — 1, lifetime effects lead to big effects in Re[E(k,io)]. For example, the curves are no 
longer regularly displaced with respect to one another. Also, the frequency range increases appreciably and the heights 
of the curves decrease. Let us comment that the numerical form of the real part of the self-energy clearly shows that 
the Kramers-Kronig relations for the self-energy are satisfied, in order to have the roots of the one-particle Green 
function on the same side of the complex plane. Analitically, it can be proved too. The only requirement coming out 
of these calculations, with our Ansatz, is that a(k) > |5lJ (see Appendix A). As a(k) x "f(k) < 0, then j(k) < 0. 
This considerations we have checked in further numerical calculations with lattice sizes of 64 x 64. Due to these new 
findings, we correct the results of Ref. Q since the solutions we found there must satisfy the conditions established 
here, i.e., a(k) > and 7(fc) < 0. 

Figs. 4a, 4b and 4c show the imaginary part of the double-occupied Green function, — Im[G2(k, lo)], vs lo along the 
diagonal of the Brillouin zone for the same values of interaction as before. Let us recall that G^/c, lo) is given by 



G 2 (k,u;) ee TV (rWr);4 (0) \ ^ , (9) 



(k,u 



where a = —a. In Eq. ([)]) (k, lo) means the Fourier transform of the spatial-temporal correlation function and T T 
the usual time ordering of the operators. Using the equation of motion technique for the one-particle Green function, 
G(k,Lo), we obtain that G(k,co) and G2{k,io) are related as follows [[32| 

(uj-e^)G(k,uj) = l + UG2(k,u J ) . (10) 

We observe that there is a big peak in the interval lo G [—2, —1] which is most likely due to the peak in A{k,Lo). 
However, the right peak at (k,ui/W) = (it, tt, w 1/3) increases with interaction. At the same time, we see that the left 
frequency peaks (to < 0) start to line up for small momenta but they almost vanish for k — 7r(3/4, 3/4) and 7r(l, 1). 

X (2) {k,io) = — km Im\G 2 (k,w + i8)] , (11) 

is the spectral density for the double-occupied Green function. We see from Figs. 4 that there is negative contribution 
to this spectral density, which is due to the presence of the factor lo — in front of the one-particle spectral density, 

A(k, lo). In addition, the factor ui — is controlling the height of the peaks in x^ik, lo). For example, when 

X^ 2 \k,Lo) — is because this factor is zero. x^ 2 \k,ijo) is given by 

U X {2) (k,Lo) = (uj-e s )xA(k,Lo) , (12) 
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which is identically zero for the non-interacting electron gas, since A(k 7 to) is a Dirac delta function at the same 
argument of the quantity in front of it. So, any deviation from zero is a signature of an interacting system. Contrary 
to A(k,uj), which is always positive, \ can be negative. The only requirement is that p2| 

/+oo 
X (2) (k,u;)dLo = p , (13) 
-oo 

where p is the electron density/spin. This can be easily checked calculating the first moment or moment of zeroth 
order for the double-occupied Green function. The relation between the self-energy and the double-occupied Green 
function is the following 

t/G 2 (fc,c) = 7 S(fc '^ . (14) 

Eq. ( |l4| ) is an exact relationship and it can be used to keep control of th e ap proximations made in the self-energy 
and the double-occupied Green functions, as it has been discussed in Ref. |^3[. Needless to say that to approximate 
Y,(k,u) is equivalent to have an approximation for G2(k,ui) and viceversa. In consequence, simple approximations 
for G2(k,u) are not always equivalent to simple approximations for G(k,u>) (or S(fc,cj)) or viceversa. For example, 
a single pole ansatz (without lifetime effects) in G2{k,uS) leads to the Hubbard-I solution as it has been discussed in 
Ref. Q. To go beyond the Hubbard-I solution for G 2 (k,ui) we have to use Eqs. dlpl). 



IV. CONCLUSIONS AND FUTURE TRENDS 

We have investigated the dynamical quantities, A(k,ui), Re[Y.(k, u>)], —Im[T,(k,uj)] and —Im[G2(k,uj)], vs to along 
the diagonal of the Brillouin zone, for three values of the interaction, namely, U/W — 1/2,2/3,1. In all these 
quantities we observe that the role of correlations and lifetime effects is fundamental. For example, for values of 
U/W « 1 the one-particle spectral density becomes almost one-peak, while — Im[E(k, lo)] becomes a wider inverted 
Lorentzian. Re[E(k,ui)}, for U/W — 1 has lost all sign of regularity. A(k,u>) becomes featureless for large values of 
U/W. Our treatment of G(k, u) and G2(k, u>) is not perturbative since we impose sum rules to A(k, lu) to find E(fe, uj) 



and G 2 {k, lo) is found from the equation of motion technique (Eq. (10)). 

The choice of self-energy (Eq. (||) is an attempt to shed some light on Nolting approach to which Eq. (||) reduces 
when j(k) — fis}| . Nolting's study (when looked upon with our optics, i.e., j(k) = 0, in Ref. JlJ]) is also a non-Fermi 
liquid. Our Ansatz for E(fc, u)) is rather phenomenological, since we have not invoked any microscopic mechanism to 
postulate it (Eq. (|^)). However, we have been guided b y t he single pole structure of Nolting without lifetime effects. 
This structure has been fleshed out in a recent paper p8[ . Also, we have relied on the calculations of Kishore and 
Granato |2C| ] which represent a non-Fermi liquid approach for the self-energy. Those interested in a nice interpretation 
of non-Fermi liquid behavior of the experimental data of High-Temperature Cuprates, please see Ref. [Q. Work is 
in progress p5l to include Fermi liquid features close to the chemical potential. According to our belief, this type of 
considerations are much harder to be tackled with the procedure of Nolting, i.e., two poles in the one-particle spectral 
function, A(k,u>). For example, the self-energy proposals of Norman et al [ p6| , for the overdoped and underdoped 
regimes of the cuprate superconductors, can be numerically solved for the attractive Hubbard model J37J, for d-wave 
superconductivity, where off-diagonal Green functiond are called for. 

A single pole structure in E(fc, u>) goes beyond the Hubbard-I approximation, since the Hubbard-I approximation 
is also equivalent to choose a single pole in G2(k,u>) (without lifetime effects). This can easily be checked due to the 
exact relationship given in Eq. (J14|). However, by comparing the results for the moments without lifetime effects we 
find that a(k) « p(l — p)U 2 which proves that our choice for E(fc, w) is, at least, a second order expansion in U. This 
is in agreement with the theoretical findings of Appendix A, since a(k) > 0. Thus, a(k) is almost fc-independent. 
Similarly, we find that j(k) < is independent of k, but strongly dependent on U . In consequence, our numerical 
study proves that our Ansatz is the easiest way to include lifetime effects and to consider Fermi and/or Marginal 
Fermi liquid behavior in the original proposal of the moment approach of Nolting. 

We thank CONICIT- Venezuela (project F-139), the Brazilian Agency CNPq, FAPERGS and the Swiss National 
Science Foundation for finantial support. Interesting discussions with Prof. H. Beck, Prof. M.S. Figueira, Prof. E. 
Anda and Dr. M.H. Pedersen are fully acknowledged. In particular, Prof. Beck brought to our attention Ref. [ p6[ . 
Thanks to Maria Dolores Garcia Gonzalez for a reading of the manuscript. 
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V. APPENDIX A: POLES OF THE ONE-PARTICLE GREEN FUNCTION 

With the self-energy ansatz given by Eq. (|J), the one-particle Green function becomes 

G(k, U ) = *-n(fc)-ry(fc) ^ (lg) 

(w - e £ )(w - Q(k)) - a(k) - tj(k)(uj - eg) 

From Eq. ( |l5| ) the poles of the one-particle Green function are given by the roots of the following equation: 

z 2 -{E]: + n{k)+ij(k))z-a(k) + (n(k)+i"f(k))£]: = . (16) 

Solving Eq. (|l^) we get that the two roots are 



eg + (l(Je) + i"f{k) ± J{e i -Sl{k)-i 1 {k)y +Aa{k) 
*± = 2 (17) 

We have to find the real and imaginary parts of the two roots. For this we follow the standard procedure making 



y/x + iy = xi + iyi , (18) 

from where we get that 

I \/x 2 + y 2 + x\ ( x/x 2 +y 2 - x\ 

X '={ 2 J ; 2/1 ={ 2 J ' (19) 

Comparing Eqs. @ and (|TJ) we conclude 

x= ( £ ~-n(k)) 2 +4a(k)-j 2 (k) ; y = 2 7 (fc)(n(*) - eg . (20) 

In consequence, z± are given by 

z± = ^(fc)+eg±q :i +z(7(fc)±;/ 1 ) (2i) 
If we require that our roots be on the upper half-complex plane, we must impose that j(k) ±y\ > 0. Carrying out 



the calculations we arrive to the result that a(k) > 0, which proves the statement advanced in Section III 
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FIGURES 



Figs la, lb and lc. A(n, n, w) vs u) along the diagonal of the Brillouin zone for three different values of interaction, 
namely, U/W = 1/2, 2/3 and 1. Our system has a periodicity of 32 x 32. We arc at half filling, p = 1/2. As we work 
in two dimensions, the bandwidth is W = 8t. The wave vector along the diagonal is defined as k = 

Figs. 2a, 2b and 2c. — 7ro£(n, n, oj)] vs u along the diagonal of the Brillouin zone. Same parameters of Fig. 1. 

Figs 3a, 3b and 3c. Re[T,(n, n, u)] vs n along the diagonal of the Brillouin zone. Same parameters as previously. 

Figs. 4a, 4b and 4c. — Im[G2(n, n, u))] vs uo along the diagonal of the Brillouin zone. Same parameters as before. 
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